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Abstract 

We present the fortran code SusyBSG version 1.1, which computes the branching ratio for 
the decay B 7 in the MSSM with Minimal Flavor Violation. The computation takes 

into account all the available NLO contributions, including the complete supersymmetric QCD 
corrections to the Wilson coefficients of the magnetic and chromomagnetic operators. 
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1 Introduction 



The radiative B decays play a key role in the program of precision tests of the Standard 
Model (SM) and its extensions. The inclusive decay B ^ Xg'y is particularly well suited to 
this precision program, thanks to its low sensitivity to non-perturbative effects. The present 
experimental world average for the branching ratio of this decay, with a 1.6 GeV lower cut on 
the energy of the photon, is BR[B X, 7]cx = (3.55 ± 0.26) x 10"^ [Ij. The SM prediction for 
the branching ratio with the same cut on the photon energy is BR[i? — >■ Xs 7]th = (3.15±0.23) x 
10""^ [21 [3] and includes most of the next-to-next-to-leading order (NNLO) perturbative QCD 
contributions as well as the leading non-perturbative and electroweak effects. Both experiment 
and SM prediction have an uncertainty of about 7%. 

New Physics (NP) can in principle induce sizeable contributions to the decay B —>■ Xsj, 
hence the good agreement between the SM prediction and the experimental result puts severe 
constraints on the flavor structure of NP models. However, the theoretical accuracy of the 
predictions for BR[i? Xg 7] in extensions of the SM is not at the same level as in the SM. A 
complete next-to-leading order (NLO) calculation is available only for the Two-Higgs-Doublet 
Model (THDM) [UlS], the Left-Right symmetric model [6], and for the Minimal Supersymmetric 
Standard Model (MSSM) under the simplifying assumption of Minimal Flavor Violation (MFV), 
according to which the quark and squark mass matrices can be simultaneously diagonalized 
and the only source of flavor violation is the CKM matrix. In this scenario, the diagrams that 
include gluons and charginos were computed in refs. O E] , while the contributions involving 
gluinos were first considered in the heavy gluino limit in ref. [7], and in an effective Lagrangian 
approach in refs. P, El [ID!- After a partial two-loop calculation [TT], the full computation of 
the two-loop gluino contributions to i? ^ 7 in the MSSM with MFV was finally presented 
in ref. [12j. 

Several public computer codes that determine the MSSM mass spectrum and other SUSY 
observables (e.g. SuSpect [I3], SPheno [H], micrOMEGAs [l5], FeynHiggs [16], NMHDECAY [TTj . 
CPsuperH [TB] and Superlso [12]) contain calculations of BR[i? Xg'j] in various approxima- 
tions. However, in the present versions of all these codes the two-loop gluino contributions to 
i? — > Xs 7 are included, if at all, only in the effective Lagrangian approximation of refs. [8tl9|[T0]. 
which is valid in the limit of heavy superpartners and large tan/3. 

In this paper we present a new fortran code, SusyBSG, dedicated to the full NLO calculation 
of BR[B Xs 7] in the MSSM with MFV. The code includes the full results of ref. [12] for the 
two-loop gluino contributions to the Wilson coefficients of the magnetic and chromomagnetic 
operators relevant to the S — X^ 7 decay, and the results of refs. [3 |6] for the two-loop gluon 
contributions. It should be recalled that the weak interactions affect the squark and quark mass 
matrices in a different way, therefore their simultaneous diagonalization can be consistently 
imposed only at a scale /Xmfv, which concurs to specify the MFV model. The renormalization 
group evolution of the MSSM parameters then leads to a disalignment between the squark and 
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quark mass matrices at scales different from /Imfv Thus, optionally, the code allows for the 
inclusion of (small) additional contributions to the Wilson coefficients from one-loop diagrams 
with gluinos and down-type squarks (as well as charginos and up-type squarks) that occur 
when the MFV condition is imposed at a scale much higher than the weak scale. For the sake 
of comparison, SusyBSG can also provide evaluations of BR[i? — »• Xg 7] in the SM and in the 
THDM, as well as in the MSSM with two-loop gluino contributions computed in the effective 
Lagrangian approximation. 

In SusyBSG the relation between the Wilson coefficients and the B —>■ Xs'j branching ratio is 
computed at NLO in perturbative QCD, along the lines of ref. [20], including also the dominant 
electroweak and non-perturbative corrections. For new physics that does not induce effective 
operators other than those already present in the SM, the NNLO anomalous dimensions of 
the effective operators and their matrix elements are the same as in the SM, suggesting the 
possibility of a partial NNLO implementation. A complete NNLO calculation would require 
also the NNLO contributions to the matching conditions, which have been computed only in 
the SM. However, if it can be argued that these contributions are as small as in the SM they 
can safely be neglected. This has been done for the case of the type II THDM in ref. [2], where 
the NP contributions to the matching conditions are computed at NLO and the anomalous 
dimensions and matrix elements are computed at NNLO. In the present version of SusyBSG 
we follow a similar but simpler route: we take into account the first NNLO estimate of ref. j2j 
by modifying the approach of ref. [2Dj in a way that approximately reproduces the results of 
a partial NNLO implementation. In fact, any attempt at a partial NNLO implementation 
has limitations in the MSSM, where higher-order QCD and electroweak contributions to the 
matching conditions can be sizeable, especially at large tan (3. The theoretical accuracy of our 
code therefore tends to be poorer in the MSSM than in the SM. 

This manual is structured as follows: in section 2 we briefly summarize the two-loop results 
implemented in SusyBSG, focusing on the information necessary to a correct interpretation of the 
input parameters (we refer the readers to ref. [12j for the technical details on the calculation). In 
section 3 we describe the structure of SusyBSG, focusing in particular on the input and output 
parameters of the two main subroutines that make up the program. In section 4 we briefly 
detail our default choices for the input parameters and discuss the theoretical uncertainty. In 
the appendices we provide additional details on the various corrections implemented in SusyBSG. 
In the appendix A we present the one-loop gluino contributions to the matching conditions for 
the Wilson coefficients Ci^^ and C2^\ In the appendix B we provide formulae for the one-loop 
gluino and chargino contributions to the Wilson coefficients in the presence of flavor mixing 
in the squark sector. In the appendix C we summarize our treatment of the tan /3-enhanced 
contributions to the Wilson coefficients. Finally, in the appendix D we summarize the NLO 
computation of the relation between Wilson coefficients and branching ratio for B —>■ Xsj. 

The latest version of SusyBSG can be downloaded from the Web page 

http : / / cern . ch/ slavich/ susybsg/home . html 
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2 Radiative B decays in the MSSM with MFV 



In this section we summarize the calculation of the weak-scale matching conditions for the 
AB = 1 effective Hamiltonian in the MSSM with Minimal Flavor Violation, as implemented in 
SusyBSG. The NLO relation between the Wilson coefficients at the weak scale and the branching 
ratio is summarized for completeness in the appendix D. 

The AB = 1 effective Hamiltonian at the matching scale (of the order of the weak scale) 
is given by 

AG J? 



where Gp is the Fermi constant and Vj^^^, V^^^^ are elements of the CKM matrix. The oper- 
ators relevant to our calculation are 
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When the QCD corrections are considered, the Wilson coefficients of the operators Qi can 
be organized in the following way 
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where the various leading order (LO) contributions are classified according to whether the 
corresponding diagrams contain only SM fields, a physical charged Higgs boson and an up-type 
quark, or a chargino and an up-type squark. The expressions for C^^^^^ and Cf''^ can be 
found, e.g., in ref. while those for (^^susy found, e.g., in eq. (4) of ref. f7J. Note that 

at LO only the coefficients of the magnetic and chromo-magnetic operators Qj and receive 
contributions from non-SM fields. One-loop neutralino- and gluino-exchange diagrams should 
be neglected under the MFV assumption. 
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The NLO coefficients c^^"^^^ and Cf^'*'^ contain the gluonic corrections to the SM and 
charged Higgs contributions, respectively, and can be found for instance in ref . |i4j . Concerning 
the NLO supersymmetric contributions C^^^ susy^ ^j^^ chargino-gluon contributions can be found 
in refs. [71 [6] , and a full two-loop computation of the gluino contributions to Cy^g ^^^^ in the 
MSSM with minimal flavor violation was more recently presented in ref. [12]. Together with 
the one-loop gluino contributions to C[^l ^^^^ that we present in the appendix A, the results of 
ref. [12] provide us with a complete NLO computation of the supersymmetric QCD contributions 
to BR [5 Xs 7] in the MFV scenario. 

The study of NLO contributions in the MFV scenario is complicated by the fact that the 
simultaneous diagonalization of squark and quark mass matrices is not preserved by radiative 
corrections. As a result, even when the renormalized mixing matrices for quarks and squarks 
are assumed to be flavor-diagonal, flavor-changing counterterms of electroweak origin have to 
be taken into account in the two-loop calculation. The technical issues related to the renor- 
malization of flavor mixing are discussed in ref. [12] , to which we refer the interested reader. In 
the following we summarize the information necessary to a correct interpretation of the input 
parameters required by SusyBSG. 

In the MFV framework the mass matrices for both the up-type and down-type quarks and 
squarks are assumed to be simultaneously flavor-diagonal at some renormalization scale //mfv 
When computing the matching conditions for the Wilson coefficients we neglect the masses of 
the first- and second-generation quarks, as well as the left-right mixing in the mass matrices 
of the corresponding squarks. Thus, the soft SUSY-breaking terms that enter the squark mass 
matrices and are relevant to our calculation are: the masses for the SU(2) squark doublets, 
mg., where i is a generation index; the masses for the third-generation singlets, m-p and m^; 
the trilinear interaction terms for the third-generation squarks. At and Ah. We recall that 
in the so-called super-CKM basis, where the matrices of Yukawa couplings are diagonal and 
the squarks are rotated parallel to the quarks, the 3x3 mass matrices for the up-type and 
down-type left squarks are related by {M^)ll = V^^^ {M^)ll ^cKMt_ Therefore, the two mass 
matrices can be both flavor- diagonal only if they are flavor-degenerate. This means that the 
MFV scenario can be consistently implemented only if we choose a common mass parameter 
for the three generations of SU(2) squark doublets, i.e. mq. = mq. 

Since we are focusing on the QCD corrections at two loops, it is necessary to specify a 
renormalization scheme for the input parameters mg, m-p and At, which determine the up- type 
squark masses and mixing entering cf^^^^^ and are subject to 0{as) radiative corrections. We 
consider two options: the ffist is to assume that they are Lagrangian parameters expressed in 
a minimal subtraction scheme such as DR, at some renormalization scale /^susy of the order 
of the superparticle masses. The second is to adopt an on-shell (OS) definition: mg, rriT and 
At can be interpreted as the unphysical parameters that enter the tree-level stop mass matrix 
obtained by rotating the diagonal matrix of the physical stop masses by a suitably defined 
physical mixing angle Oi- In this case we adopt an OS definition (i.e. we use the physical mass) 
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also for the top quark mass that enters the tree-level stop mass matrix. 

The other MSSM parameters relevant to the calculation, for which we need not specify a 
renormalization prescription, are: the ratio of Higgs vacuum expectation values tan/3 = f2/fi; 
the charged Higgs boson mass m/^±; the gluino mass mf, the SU(2) gaugino mass parameter 
M2; the higgsino mass parameter /i. Our conventions for the signs of the various Lagrangian 
parameters are as specified by the SUSY Les Houches Accord (SLHA) [2T1 i.e.: the top 
and bottom quark masses are nit = htV2/V^ and = h^V\l \f2, respectively; the left-right 
mixing terms in the stop and sbottom mass matrices are mt{At — n cot j3) and mb{Ai, — tan (3) , 
respectively; the (2,2) entry of the chargino mass matrix is /x. 

SusyBSG has an option to read the parameters of the MSSM Lagrangian from a SLHA 
spectrum file. In that case, the parameters are meant to be expressed in the DR renormal- 
ization scheme. When a SLHA spectrum file is available the program reads from it also the 
DR-renormalized top quark mass computed at the scale /Xsusy, which is necessary for the com- 
putation of the running stop masses. In the absence of a SLHA spectrum file the running top 
mass is computed internally by the program, taking into account only the corrections controlled 
by the strong gauge coupling. 

In a two-loop calculation the interpretation of the MFV requirement itself depends on the 
way we renormalize the flavor mixing, i.e. the way we flx the counterterms that cancel the 
divergences of the antihermitian parts of the quark and squark wave-function-renormalization 
(WFR) matrices. In particular, if we perform a minimal subtraction we are imposing the MFV 
condition on the DR-renormalized parameters of the Lagrangian evaluated at the scale /Xmfv 
In this scheme C7^g^'^^^ contains logarithms of the ratio Ms/ ^mfv, where Mg represents the 
mass of the superparticles entering the loops. An alternative option consists in subtracting also 
the finite part of the antihermitian WFR: this results in a conventional (and gauge-dependent) 
on-shell renormalization scheme [23], in which Cy^g^^^^ is independent of /iMFv 

From the discussion above it is clear that, in models where the MFV condition is imposed 
at a scale much larger than the superparticle masses (such as, e.g., supergravity models where 
one identifies /^mfv with the GUT scale), the Wilson coefficients computed in the minimal 
subtraction scheme contain very large logarithms of Mg/n^^y. In this case, no matter what 
renormalization scheme is chosen, the fixed-order calculation does not provide a good approxi- 
mation to the correct result. Indeed, in such models the soft SUSY-breaking mass parameters 
- which are fiavor-diagonal at the scale /Xmfv ^ niust be evolved down to Ms with the appropri- 
ate RGE, thus generating some fiavor violation in the squark mass matrices. When the squark 
mass matrices are diagonalized, the resummed logarithms of the ratio Ms/ /iMFv are absorbed in 
the couplings of the resulting squark mass eigenstates with the gluinos (and the charginos and 
neutralinos). Typically, the effects of the RGE- induced fiavor mixing are relatively small, and 
we include them only at leading ordei0, in the one-loop diagrams with gluinos and down-type 

^We follow the approach of ref. 24J. For a more complete treatment of the gluino contributions in scenarios 
with generic sources of flavor violation see ref. |25j. 
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squarks (which would vanish if the MFV condition was vahd at the low scale Mg) and in the 
one-loop diagrams with charginos and up-type squarks. The corresponding contributions to 
the Wilson coefficients are given explicitly in the appendix B. Once the logarithmic effects have 
been taken into account in this way, the genuine two-loop MFV contributions in Cy^^^'^^^ can 
be computed by setting artificially /Xmfv ~ Mg. In this case, using either the on-shell scheme or 
the minimal subtraction scheme to renormalize the flavor mixing will give basically the same 
result. 

We remark here that the rather complicated task of solving the system of RGE equations for 
the soft SUSY-breaking parameters, taking into account the full flavor structure of the MSSM, 
is not performed by SusyBSG. In order to obtain the squark masses and mixing matrices that 
enter the expressions in appendix B, starting from flavor- universal boundary conditions at the 
scale /Imfv, the users can run one of the public spectrum calculators that include a full treatment 
the flavor structure, such as the latest (still unpublished) versions of SPheno [T3] or Sof tSusy 
[SB] , and pass the results to SusyBSG by means of a SLIIA2 [22] input/output ffie. 

Finally, it is well known that in the MSSM the relation between the bottom quark mass rrib 
and the bottom Yukawa coupling hb is subject to tan /3-enhanced threshold corrections [27j. If 
the bottom Yukawa coupling entering the one-loop part of the Wilson coefficients is expressed 
in terms of the SM value of the running bottom mass, the SUSY threshold corrections induce 
counterterm contributions that, although being formally of higher order in the loop expansion, 
are enhanced by powers of tan (3 and may therefore be sizeable when tan f3 is large. As discussed 
e.g. in ref. [28], such potentially large corrections can be absorbed in the one-loop results by 
a suitable redefinition of the bottom Yukawa coupling. Other tan /3-enhanced contributions to 
the Wilson coefficients appear in the form of corrections to the Higgs-quark-quark vertices, and 
have been computed in refs. [SlISllTn] in an effective Lagrangian approach where the heavy SUSY 
particles are integrated out of the theory. Among the tan /3-enhanced two-loop contributions, 
those controlled by the strong gauge coupling are fully accounted for by the calculation of 
ref. [12]. In addition, we include in SusyBSG the tan /3-enhanced contributions controlled by 
the top and bottom Yukawa couplings, following the approach of ref. [TU] . More details on the 
treatment of the tan /3-enhanced contributions are given in the appendix C. 
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3 Structure of the program SusyBSG 



SusyBSG 1 . 1 is structured in two modules, which in principle could be used independently as 
stand-alone programs. The first module is the subroutine WilsonCoef f , which computes the 
weak-scale matching conditions for the Wilson coefficients of the Ai? = 1 effective Hamilto- 
nian in a physics model to be chosen among the SM, the Two Higgs Doublet Model and the 
MSSM. The second module is the subroutine getBR, which computes BR[i? Xg'y] at NLO, 
taking as input the SM and new-physics contributions to the Wilson coefficients evaluated by 
WilsonCoeff . The main program BSGAMMA, which can be freely modified according to the 
users' needs, sets all the relevant input parameters, calls the two subroutines and prints out 
the output. Within BSGAMMA the users are allowed to read the input parameters required by 
WilsonCoeff from a spectrum file written in the SLHA format. More specifically, the sub- 
routine readSUSY_SLHAl reads from a SLHAl [21] file only the flavor- conserving parameters 
required for the calculation of the one- and two- loop contributions appearing in eq. flTU]) . The 
subroutine readSUSY_SLHA2, instead, reads from a SLHA2 [22] file all the parameters of the 
MSSM Lagrangian, including the flavor-violating parameters necessary to compute the addi- 
tional one-loop contributions discussed in the appendix B. In either case the spectrum file must 
be named SLHA . in, and it must be located in the directory where the program is run. 

For certain choices of the input parameters it may happen that the masses of two particles 
are accidentally very similar to each other, or that the sum of two masses is very close to a third 
mass. In some pathologic cases this leads to numerical instabilities in the output of SusyBSG, 
even if all the formulae for the two-loop gluino contributions to B ^ Xgj are well behaved when 
the corresponding limits are taken analytically. Since those formulae are very long and depend 
in a complicated way on the values of twelve different particle masses, providing analytical 
results to cover all the problematic limits would be highly inefficient in terms of size and speed 
of the code. Therefore, we limit ourselves to looking for the occurrence of accidentally similar 
masses, and issuing a warning if necessary. We leave it to the users to check that, in those 
cases, the result for BR[i? — >■ X, 7] is not unreasonably sensitive to small variations of the 
input parameters The only excepiion is the case in which the mass of the first- or second- 
generation up-type squark m^^ is very similar to the mass of the super-strange w^s^. This 
happens inevitably if the soft SUSY-breaking parameter mg, common to both masses, is much 
larger than rriz- To avoid numerical instabilities we switch to the analytical results valid in the 
limit rriji^ = mg^ if the relative difference between the two masses is less than 1%. 

In the following we describe in detail the input and output parameters of the two main 
subroutines that make up SusyBSG. 

^Note that the numerical instabihties can be greatly reduced by compiling the program in quadruple precision. 
Details are provided in the SusyBSG Web page. 
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3.1 The subroutine WilsonCoeff 

The call to the subroutine for the matching conditions to the Wilson coefficients reads 

call WilsonCoeff (imod, scheme, muO,mususy,mumfv, 
$ msqS , mstr , msbr , msql , At , Ab , mHp , mg , M2 , mu , t anb , 
$ CISM , C7SM , C8SM , CINP , C7NP , C8NP , prob , eqmass) 

The variables in the ffist two lines of the call command are inputs, and are defined as: 

• integer imod: allows the users to choose the particle content of the theory and the 
approximation used in the computation of the two-loop contribution to the Wilson coef- 
ficients. The values 0-4 correspond to 

0: Standard Model fields only. The new-physics contributions are set to zero; 

1: Two Higgs Doublet Model. The only additional contributions come from diagrams 
with a charged Higgs if^; 

2: MSSM with the effective Lagrangian approach. At NLO, only the tan /3-enhanced 
SUSY contributions are included in the Wilson coefficients as in ref. |10j ; 

3: MSSM with full NLO QCD contributions. Includes the two-loop calculation of the 
gluino contributions presented in ref. [12] , and the results of ref. [10] for the remaining 
(i.e., non-QCD) tan /3-enhanced contributions. 

4: MSSM with full NLO QCD contributions (as for imod = 3), with in addition the 
contributions of the one-loop diagrams with gluinos and down-type squarks and 
chargino and up-type squarks given in the appendix B. 

• logical scheme (2): contains two logical switches (relevant only for imod = 2,3,4) 
that allow the users to specify the renormalization conditions in the squark sector. In 
particular: 

1: choice of renormalization scheme for the input squark mass parameters 

(.true. = OS, .false. = DR); 
2: choice of renormalization scheme for the MFV condition 

(.true. = OS, .false. = DR). 

If the MSSM Lagrangian parameters are read in from a SLHA file both entries of scheme 
should be set to .false. 

• real*8 muO: renormalization scale fiQ (of the order of the weak scale) at which the 
matching of the Wilson coefficients is performed. 
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• real*8 mususy: renormalization scale /Xgusv at which the input squark mass parameters 
are given (relevant only for schemed) = .false.). 

• real*8 mumf v: renormalization scale /Imfv which the MFV condition is imposed (rel- 
evant only for scheme (2) = .false.). As explained in section 2, the results of the 
program are not reliable if /Imfv is set to be much larger than //susy- 

• real*8 msq3,mstr ,msbr: third-generation squark mass parameters ttiq^, uit and tub- 

• real*8 msql: mass parameter mq for the first- and second-generation squark doublets. 
Note that the MFV condition is only consistent with uiq — uiq^. 

• real*8 At,Ab: third-generation Higgs-squark-squark interaction terms At and Af,. 

• real*8 mHp: mass of the charged Higgs boson m^±. 

• real*8 mg: gluino mass rrig. 

• real*8 M2: SU(2) gaugino mass parameter (enters the chargino masses). 

• real*8 mu: Higgs-mixing superpotential parameter /i. 

• real*8 tanb: ratio of Higgs vacuum expectation values tan/5. 

of course, all of the SUSY input parameters are relevant only for imod = 2,3,4, while mHp and 
tanb are relevant for imod = 1,2,3,4. The option imod = 4 requires the presence of a SLHA2 
spectrum file, from which the subroutine readSUSY_SLHA2, to be called before WilsonCoef f , 
reads the SUSY input parameters. For consistency, when imod = 4 the entries of scheme should 
be both .false., and /Xmfv should be set close to /isusY- 

The variables in the third line of the call command are outputs, and are defined as: 

• real*8 CISM(6): vector containing the Standard Model contributions to Cf^^(//o) for 
1 < i < 6 (they are actually different from zero only for i = 1, 4). 

• real*8 C7SM(2) ,C8SM(2): vectors whose two elements are the Standard Model contri- 
butions to C7°8(/xo) and c!^^l{iJ,o). 

• real*8 CINP(6): vector containing the new-physics contributions to C^^\ij,q) for 1 < 
i < 6 (in the MSSM they are different from zero only for i — 1, 2, 4). 

• real*8 C7NP (2) , C8NP (2) : vectors whose two elements are the new- physics contributions 
to Cj^sVo) and C^^liiio). 
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• logical prob: problem flag, set to .true, if the relative difference between two (or 
more) particle masses is less than 1%. The case m^^ ^ - which is taken care of 
internally - is not flagged, nor are some other cases that never lead to instabilities. 

• logical eqmass(12): vector that identifies the masses that are too close to each other. 
The ordering is {nit, rriw, rrig, mH±, m^^, mi^, m^^, m^^, m^i^, m^^, m^+, fn^+), and the 
entries of eqmass corresponding to the masses that are too close are set to .true. 

In addition to the input variables passed in the call to the subroutine, the users must set 
(or read from a SLHA spectrum file) a number of Standard Model parameters in the common 
block SMINPUTS 

common/ SMINPUTS/ mz , mw , mtpole , mbmb , hsm , asmz , azinv 

The variables in SMINPUTS are defined as 

• real*8 mz,mw: physical masses and for the SM gauge bosons. 

• real*8 mtpole: physical top-quark mass ni\"^'^ . 

• real*8 mbmb: running bottom mass m^^(mfe) in the MS scheme at the scale mfe. 

• real*8 hsm: mass mn of the SM-like Higgs boson (used only for the small electroweak 
corrections of ref. [21], see the appendix D). 

• real*8 asmz, azinv: strong and inverse electromagnetic couplings as(mz) and l/a(m2) 
in the MS scheme at the scale ruz- 

3.2 The subroutine getBR 

This subroutine computes BR[i? — > Xs'i\ at NLO using the results of refs. |20l [29]. For the 
users' convenience the computation of the branching ratio will be summarized in the appendix 
D. The call to the subroutine reads 

call getBR (muw , mut , mub , muc , EO , CINP , C7NP , C8NP , BR) 

The variable BR is the output value for BR[i? — > Xs7]. All the other variables in the call 
are inputs, and are defined as 

• real*8 muw: weak scale ~ f^w at which the light quark contributions to the Wilson 
coefficients are computed. 

• real*8 mut: weak scale /it ~ at which the top quark contributions to the Wilson 
coefficients are computed. The new-physics contributions (see below) should also be 
computed by the subroutine WilsonCoef f at a scale /io = /ij. 
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• real*8 mub: low scale /i^ ~ at which the branching ratio is computed. 

• real*8 muc: low scale /ic at which the charm quark mass m^^(/ic) entering the SM 
contributions to the Wilson coefficients is computed. 

• real* 8 EO: minimum photon energy Eq. 

• real*8 CINP(6): vector containing the new-physics contributions to Cl^\^t) for 1 < ^ 
6 (they are produced in output by the subroutine WilsonCoef f ). 

• real*8 C7NP (2) , C8NP (2) : vectors whose two elements are the new-physics contributions 
to C^},{lJ't) and C'^^ilJ^t) (they are produced in output by the subroutine WilsonCoef f). 

In addition to the input variables passed in the call, the subroutine getBR requires the 
Standard Model input parameters contained in the common block SMINPUTS (see section [XT]) , 
and a set of parameters contained in the common block BRINPUTS 

common/BRINPUTS/ aOinv , mcmc , rbs , hlam , ccsl , bsl , lambda , A , rhobar , etabar 

The variables in BRINPUTS are defined as 

• real*8 aOinv: inverse of the fine-structure constant, l/aem- 

• real*8 mcmc: running charm mass m^^(mc) the MS scheme at the scale ttLc. 

• real*8 rbs: ratio rrib/ms entering the gluon bremsstrahlung contribution. 

• real*8 hlam: HQET parameter A2 entering the non-perturbative contribution. 

• real*8 ccsl: non-perturbative semileptonic phase-space factor C. 

• real*8 bsl: semileptonic branching ratio BR[i? ^ Xc e z/]. 

• real*8 lambda, A, rhobar , etabar: Wolfenstein parameters X , A , p and fj for the CKM 
matrix. 

More information on these quantities can be found in the appendix D. Note that the param- 
eters entering BRINPUTS are not read in from a SLHA spectrum file (even when one is present) 
and must be explicitly set by the users before calling getBR. 
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3.3 File structure of SusyBSG 



The latest version of SusyBSG can be downloaded from the program's Web page in the form of a 
file named SusyBSG_vv.tar .gz, where vv stands for the version number. When uncompressed 
and unpacked, the program consists of the following files: 

• BSGAMMA.f: contains the main program BSGAMMA, to be modified (or replaced) by the 
users. It just sets the relevant inputs, calls the two subroutines WilsonCoef f and getBR 
and prints out the output. 

• DGStwoloop . f : contains the subroutine DGSbsgamma, which computes the two-loop gluino 
contributions to the Wilson coefficients, based on the results of ref . [T2j . Upon compilation 
it includes all the files located in the directory source. 

• f unctions. f: contains the definitions of various functions entering the one-loop and 
two-loop contributions to the Wilson coefficients. 

• getBR. f: contains the subroutine getBR, based mostly on the results of ref. [20], which 
computes the branching ratio for the process i? — > Xj 7 taking as input the new-physics 
contributions to the Wilson coefficients. 

• Makefile: by default it uses the compiler g77 to produce an executable named SusyBSG. 
Further instructions for the compilation can be found in the program's Web page. 

• slha2io.f: contains two routines that read the input parameters (with the exception of 
those in the common block BRINPUTS) from a spectrum file in the SLHA format [2TI |22]. 
readSUSY_SLHAl reads only the parameters necessary to the calculation of the MFV con- 
tributions, while readSUSY_SLHA2 reads also the parameters necessary to the calculation 
of the one-loop diagrams that involve fiavor- violating gluino-quark-squark or chargino- 
quark-squark vertices. 

• SLHA. in: an example of SLHAl spectrum file that can be used as input by SusyBSG. The 
file was produced with Sof tSusy, with input parameters corresponding to the so-called 
SPSla' [30] point. 

• source: the files in this directory contain the explicit expressions of the two-loop gluino 
contributions to the Wilson coefficients, to be included in DGStwoloop. f upon compila- 
tion. The files were automatically converted to fortran from Mathematica format using 
FormCalc [31]. 

• WilsonCoef f . f : contains the subroutine WilsonCoef f, which computes the one-loop and 
two-loop new-physics contributions to the Wilson coefficients (calling DGSbsgamma for the 
two- loop gluino contributions). 
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4 Default values of the SM input parameters 

In a phenomenological study of BR[i? Xg 7] in the MSSM with MFV, the SM input param- 
eters entering the common blocks SMINPUTS and BRINPUTS will presumably be fixed once and 
for all at the beginning of the calculation. For the users' convenience we list below the default 
values of these parameters in the version 1 . 1 of SusyBSG. In most cases we use the same values 
as in ref. [3], where the users can look for the corresponding references. 

The default values of the SM parameters entering the common block SMINPUTS are: 



mz 




ruz 




91.1876 


GeV 


mw 








80.403 


GeV 


mtpole 




pole 




170.9 


GeV 


mbls 








4.2 


GeV 


hsm 








115 


GeV 


asmz 




as{mz) 




0.1189 




azinv 




l/a{mz) 




127.918 





The default values of the B-physics parameters entering the common block BRINPUTS are: 



aOinv = 1/aem = 137.036 

mcmc = mp(me) = 1.224 GeV 

rbs = mh/ms — 50 

hlam = A2 = 0.12 GeV^ 

ccsl = C = 0.580 

bsl = BR[5 ^ Xc e z/]exp = 0.1061 

lambda = A = 0.2272 

k = A = 0.818 

rhobar = p = 0.221 

etabar = V = 0.340 



The choice of the four independent renormalization scales that appear in the calculation 
of the branching ratio and are required as input by the subroutine getBR deserves a separate 
discussion. We recall that these scales are: the weak scales pt and at which we compute 
the contributions to the matching conditions for the Wilson coefficients coming from top and 
charm quarks, respectively; the low scale /ib at which we compute the matrix elements for the 
6 — >■ S7 decay; the scale /Xc at which we express the charm mass entering the matrix elements. 
We adopt the default values 

A^f = A^iv = 2 rrinr , /i^ = 2.5 GeV , /ic = 1 GeV. (11) 

Using these values for the scales, and setting nif^'^ = 171. 4 GeV, m^^{mb) =4.15 GeV and 
l/a{mz) = 128.940 in order to reproduce the inputs of ref. [2j, we obtain a SM prediction 
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for BR[i? —>■ Xg'j] of 3.15 x 10~^, in agreement with the NNLO result of ref. [2]. Very good 
agreement, within up to 2%, is also found at various values of tan/? and mH± with the partial 
NNLO implementation of the THDM from refs. [21 [32]. While we plan to include the known 
NNLO corrections in a future version of SusyBSG, the present implementation provides an 
excellent starting point. 

It should be noted, however, that /ic in eq. ([TT]) is adjusted to a very low value in order 
to mimic the NNLO contributions that are not present in our calculation. Therefore, in this 
case, the variation of the renormalization scales should not be used to estimate the intrinsic 
uncertainty of our calculation. Indeed, the result of the NLO calculation depends quite sharply 
on /ic around the value that reproduces the NNLO result. For example, using fic = m^^(mc) = 
1.224 GeV, with the other scales fixed as in eq. ([IID, results in BR[B ^ X,7] = 3.29 x 10"^ 

Concerning the theoretical uncertainty of our prediction, we recall that in the SM analysis 
of refs. [21 [3] the error is dominated by a 5% uncertainty due to unknown O {a sAqcd/ mi,) non- 
perturbative contribution to the matrix elements. Additional ~ 4% intrinsic uncertainty stems 
from the perturbative part of the calculation and from the estimate of missing NNLO contri- 
butions. The parametric error in the SM is only about 3%. All these errors are present in the 
MSSM calculation as well. Therefore, we recommend using at least a 7% intrinsic uncertainty, 
throughout the parameter space, to be added in quadrature with the parametric uncertainty. 
Barring the case of subtle cancellations between various supersymmetric contributions, this 
appears to be a realistic guesstimate of our theory uncertainty. 

To conclude, we quote the results of SusyBSG 1 . 1 in a characteristic point of the MSSM 
parameter space, the so called SPSla' point [30]. We used SoftSusy to generate a SLHAl 
spectrum file that provides the subroutine WilsonCoef f with the DR-renormalized SUSY pa- 
rameters evaluated at a scale /isusY equal to the geometric average of the stop masses. We set 
A^MFv to be equal to ^svsy, thus enforcing the MFV condition at the weak scale. We also set 
the SM input parameters to the default values listed above (with the exception of rriw and rrih, 
which are taken from the output of SoftSusy), and the renormalization scales appearing in 
the calculation of the branching ratio to the values given in eq. (ITTl) . For the SM prediction 
(obtained with imod = 0) we find that BR[B ^ X,7] = 3.16 x 10"^. For the THDM pre- 
diction (imod = 1) we find 3.91 x 10~^. For the MSSM prediction in the effective Lagrangian 
approximation (imod = 2) we find 2.42 x 10~^. Finally, for the MSSM prediction based on the 
full NLO QCD calculation (imod = 3) we find 2.60 x 10"^. 
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Appendix A: gluino contributions to Ci^^ and C2^^ 



A complete NLO computation of the supersymmetric QCD contributions to BR[i? X, 7] 
requires the one- loop contributions of 0{as) to the Wilson coefficients of the four-fermion 
operators Qi-e- All the contributions from Standard Model particles, as well as those arising 
from loops involving a charged Higgs boson and a quark or a chargino and a squark, can be 
found in the literature [H [3 E]. There are however additional contributions to the Wilson 
coefficients Cj^"* and C2^^ that originate from loops involving gluinos, see figs. [1] and [2j These 
contributions were neglected in earlier computations performed under the assumption of heavy 
gluinos, but they must be taken into account in a general computation. Assuming MFV, and 
neglecting rric and rUs as well as the left-right mixing in the charm- and strange-squark sectors, 
the SUSY contributions to C[^^ and 6*2^^ at the matching scale /io read 



+ 



2m| 



a 



(1) SUSY 



(Al) 
(A2) 



where: summation over the repeated indices is understood; is the mixing matrix for the 
sbottoms, defined as (&i,&2)^ = ^l-,^r)^ \ Cp = 4/3 is a color factor; the quark-squark- 
chargino couplings are defined as 



rriyy COS P 



U,2Xt-V2UaX', LI = -V2Ua, L', = -V2V* 



(A3) 



where U and V are the unitary matrices that diagonalize the chargino mass matrix according 
to U Ai^V'^ = diag (m^^^, m^j). Finally, xp = rn^p/m^ for a generic particle P, and the loop 
integrals Jj are defined as 



/i(xi,a;2,a;3) 
h{xi,X2,x^) 



Xi Inxi 



(1 - Xi){xi - X2){xi - X3) 

—x\ Inxi 

(1 - Xi){xi - X2){xi - X3) 



+ (xi X2) + {xi 
+ {Xi ^ X2) + {xi ^ X3 



(A4) 
(A5) 



l3(xi,X2] 



Xi + X2 — 2 xi X2 Xi (x? — 2 X2 + Xi X2) In xi 

H 7 7TT77 ^ + [Xl 



2(xi-l)(x2-l) 



{xi - iy{xi - X2) 



X2 



(A6) 
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Figure 1: Box diagrams contributing to the Wilson coefficient 




c s be 



Figure 2: Vertex diagrams contributing to the Wilson coefficient C2 ■ Self-energy insertions 
with a squark-gluino loop on the external quark legs are also taken into account. 
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Appendix B: contributions from squark flavor mixing 



In this appendix we summarize the well-known results for the gluino-squark and chargino- 
squark contributions to the Wilson coefficients in presence of flavor mixing in the squark sector. 
In the super-CKM basis, where the squark fields are rotated parallel to their fermionic partners, 
the mass eigenstates are linear combinations of flavor eigenstates, according to 



( u, \ 
\u, J 



UL J- UR 



cl 

h 

cr 
V J 



6/4 
\de J 



r 



DL 



DR 



( d,\ 

SL 
h 

dR 
sr 

V J 



(Bl) 



where Tql and Tqr (with q = u,d) are 6x3 mixing matrices. The contributions to the Wilson 
coefficients from one-loop diagrams with gluinos and down-type squarks, 6*7*8^, have to be 
added to the ffist hne in eq. ffTUl) . They read: 



(0)5 



9s ml 



9^ ml Vt 



CKM *T/CKM 
g • ts ^tb 



E 



{^^2 (r 



2 



-7 i^*DL)i2 (rZ)ij)i3 97,S 

IlLlj 




(B2) 



where gs and g are the strong and SU(2) gauge couplings, respectively, and the functions /r^g 
and (77^8 are defined as: 

2(2 + 5x-x^) Ax 4(1 +x) 8x 



frix) 



fsix) 



27(x- 1)3 
'11 - 40a; - 19x2) 



x(l — 9x) 

: 7T In X , 

6(a;-l)4 



97[X) 



9{x 



- 13a; 



+ 



9(x- 
x{\ 



a; 



+ 



— Inx, 



9a;) 



(B3) 



36(x-l)3 6(a;-l)4 ""^"'^ 3(x - 1)^ ' 3(x - 1)^ ^"""^ ' 

When the information about the mixing in the squark sector is available we replace the 
char gino- squark contributions to C^^g^'^^"^ appearing in eq. (ITU]) with 



(0) SUSY 

7,8 



T/CKM *T/CKM 
T/CKM *T/CKM 



E 



3^ 



' 2 



-^7,8 



+ Li* R-j 



Xi 



' 1 



(B5) 

where: the indices i,J run over the up-type quark flavors; the functions F^^g and Fyg can be 
read, e.g., in eqs. (2.4) and (3.2)-(3.3) of ref. [Sj; the chargino-quark-squark couplings Ljj and 
read 



mt, 



(r 



where the matrices U and V are defined after eq. ( lA3p . We checked that eqs. ( 1B2I) and (jBSj) 
agree with the corresponding results in ref. [21]. 



V2V:,{Tul),I: 



1 



UL)jI 



(B6) 
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Appendix C: tan /3-enhanced contributions 



In this appendix we describe in detail the treatment of the tan /3-enhanced contributions to the 
Wilson coefficients relevant tob 57. The contributions arising from the threshold corrections 
to the relation between the bottom mass and the bottom Yukawa coupling [27] can be absorbed 
in the one-loop part of the Wilson coefficients by a suitable redefinition of the Yukawa coupling 
hi,. Indicating by — Cf, the tan /5-enhanced part of the SUSY contribution to the bottom quark 
self-energy and by —61, the rest of the SUSY contribution, so that the relation between the 
running bottom masses in the SM and in the MSSM is 

^ ^MSSM^^ ^ ^^^^^^ ^ ^ (CI) 

we multiply the contributions of the one-loop diagrams that involve a bottom Yukawa coupling 
by a factor k, defined as 

' - (C2) 



1 + ebtan/3 

This amounts to defining hb in terms of the MSSM running bottom mass in the one-loop 
contributions to the Wilson coefficients, with the result that the counterterm contributions at 
higher orders do not contain terms enhanced by powers of tan f3. For consistency we multiply 
by K, also the contributions of the two-loop diagrams that involve a bottom Yukawa coupling, 
but we neglect there the small effect of Sb- The 0{as) contributions to and 6b read 

(s) _ «s(/io) 2/i I ^ . 




(C4) 



where H2{x,y) is defined in eq. (2.8) of ref. [8], and Bi^p"^ , ml, ml, fi"^) is a Passarino-Veltman 
function defined as in ref. |34j . 

In addition to the tan /3-enhanced contributions that arise from the threshold corrections to 
the bottom mass, there are tan /3-enhanced contributions arising from corrections to the Higgs- 
quark-quark vertices [8], [9l [10]. We stress that all the contributions of this kind controlled by 
the strong gauge coupling are fully accounted for by the exphcit two-loop calculation of ref. [T2] . 

The tan /3-enhanced contributions to the Wilson coefficients controlled by the top and bot- 
tom Yukawa couplings are implemented in SusyBSG using the effective Lagrangian approach of 
refs. [HI, [9l [To] . They include a contribution of 0{at) to the bottom quark self-energy 




(C5) 
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such that we must use e?, = e^'''' + e[*^ when computing the factor k in eq. (lC2p . In addition, there 
is a contribution to the Wilson coefficients arising from a correction of 0{at) to the effective 
G^-t-h vertex: 



('2) 

where the loop functions FjJ(x) are defined e.g. in eq. (2.4) of ref. [8j. In the limit of heavy 
superparticles the contribution in eq. (1C6P cancels an analogous term originating from the 
redefinition of the bottom Yukawa couphng, thus ensuring the decoupling of new-physics effects 
from the SM contribution. 

There are also contributions to the Wilson coefficients involving the bottom Yukawa coupling 
and higher powers of tan/5: 

(l + e,tan/?)(l + 6l^Uan/3) ^ *^ ^ 

tan^/? ml ^^^^ 



(1 + efe tan/3)2(l + e['^ tan/3) 36 mj,^ 
where the 0{ab) vertex correction ef ■* is defined as 



Being suppressed by an additional loop factor and by ml/m'j^±, respectively, the two terms in 
eq. flC7p tend to be numerically small, but, as stressed in ref. |10], they can become relevant in 
special cases where the leading 0{as) and 0{at) effects cancel each other. 

When implementing in SusyBSG the tan /5-enhanced contributions controlled by the top and 
bottom Yukawa coupling, eqs. fICSp - flCSp . we have neglected the mixing between superpartners, 
approximating the squark masses with the corresponding soft SUSY-breaking terms and the 
masses of the higgsino components of charginos and neutralinos with the superpotential pa- 
rameter /X. Indeed, the effective Lagrangian approach used to derive these results relies on the 
assumption that the superpartners are much heavier than the weak scale and can be integrated 
out of the theory, leaving behind only non-decoupling corrections to the Higgs-quark-quark 
vertices. In this case the effect of the mixing, which is due to the breaking of electroweak sym- 
metry, can reasonably be expected to be small, and anyway is of the same order of magnitude 
as other effects that are neglected under the effective Lagrangian approximatiorlfl. 



^Formulae for the tan /^-enhanced contributions to the Wilson coefScients that take into account the effect 
of superpartner mixing have been presented in refs. [8l 115 1 135 } I36j. 
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Appendix D: NLO determination of BR[S — > ^sj] 



In this appendix we summarize the NLO computation of the relation between the branching 
ratio for B ^ Xs'~f and the Wilson coefficients computed at the electroweak scale, as imple- 
mented in SusyBSG. We rely mainly on the results of ref. pOj, but we also take into account the 
matrix elements for the four-fermion operators computed in ref. [37] (resulting in an updated 
table of "magic numbers") and the electroweak matching conditions computed in ref. |29j . 

The B ^ Xs'y branching ratio with a photon energy cutoff Eq in the 5-meson rest frame 
can be related to the experimentally measured rate BR[i? XcCujex^p as 



BR[B ^ Xs-f]E,>Eo 
BR[B ^ Xc e z/] 



cxp 



T/CKM *T/CKM 



T/CKM 
^Cb 



6 ttp 



TT C 



Kf + B{Eo) + N{Eo] 



(Dl) 



where C is the "non-perturbative semileptonic phase-space factor" defined in the appendix C 
of ref. [20]. The quantities B{Eo) and N{Eq) represent the gluon-bremsstrahlung and non- 
perturbative contributions, respectively, and will be discussed below. The factor K contains 
the contributions to the b s'j amplitude and is dominated by the effective Wilson coefficient 
for the magnetic operator at the low scale. It can be written as 



K = K, + r{fxt) [Kt + i^Np] + e. 



(D2) 



Following ref. [20] we separate K into the light quark contribution K^.^ the top-quark contribu- 
tion Kt, the new-physics contribution K^^ and the electroweak correction ^ew In the top and 
NP contributions we keep the bottom Yukawa coupling renormalized at the matching scale /it 
by introducing the quantity r(/it), defined as 



m 



MS 



m 



IS 



m 



MS 



X 



a 



Wb ) 



16 , m 
— + 81n- 

. 3 /ife 



IS\2 




(D3) 



where the "IS mass" ml^ is defined as half of the perturbative contribution to the T-mass. We 
extract the value of ml^ from the input value of m^^(mf,) using eq. (168) of ref. [38] • Unlike 
ref. [20], we evaluate r(/xt) strictly at NLO. Indeed, in eq. (1D3P the ratio r{fit) is expressed 
as the product of a term that accounts for the evolution of the running bottom mass from fit 
to Hb and a term (in the square brackets) that accounts for the shift between m^^(/ifc) and 
ml^. The NLO relation between the values of a running quark mass evaluated at two different 
renormalization scales and fiQ reads: 



m 



MS 



m 



MS 



asifi) 


TO _ 

2/3o 


_as(/io). 





1 + 



«s(/io) 7o 
47r 2(3o 




as(/i) 



(D4) 



where Pq = 11 



■ n 



f ' 



Pi = 102 — y rij , 7o = 8 and 7i = ^ — y % ('^z being the number 



of active quark flavors). The strong gauge coupling as(/i) is extracted at NLO from ag 
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according to 



as{mz) ( ^ (3ias{'mz)\nv'' 

= 1 - T, 1 

V \ /5o 47r V ^ 



v = l + f3o— In 



(D5) 



The light-quark contribution to K reads 



23 i| 8 / 11 i| 



16\ ° 

' k=l 



w 



+ 



+ 



EC 



An 

■^CKM *-^CKM 



k 

^" fc=i 



46 



15 



Ofc c?fc In h r^w In 



ub 



T/CKM *T/CKM 



An 



+ dk + dkr]yy + a{z) + rf^ 6(2;) 



(D6) 



where: /i^ is the matching scale (not necessarily equal to /ij; rjyy = as{fiw)/as{fib)] the "magic 
numbers" hk are given in eq. (2.3) of ref. [20]; the "magic numbers" ak,dk,dk,dj^ and d^ are 
given in table 2 of ref. [37]; the functions a(z) and b{z) are given in eqs. (D.l) and (D.2) of 
ref. [20]. The variable z is defined as {m^^ [fic] / fnl^^ , where fic is an arbitrary renormalization 
scale. 

After factoring out r(/it), the top-quark contribution to K reads 



l--a,(mf)2 



^ (-7(0) SM 



top iPt 
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An 



+r), 



["■k + - 

I fc=i 



„ /12523 7411 2 2 4 , ml^^ 
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'^7 topl/^tj 
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50092 1110842 16 , 32 , ml^ 
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3 ^« ^ 1^^5707^ 14283 d"" 



16 J m^"^' 
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8 top (A'-t 



(0) SM, 



(D7) 



where: rjt = as{fit)/o:s{fib)] the "magic numbers" are given in table 2 of ref. ^37j; the top 
quark contributions to the matching conditions for the Wilson coefficients at the scale fit can 
be expressed in terms of the functions Eq,Aq-^ and Fq^ appearing in eqs. (6), (10)-(12) and 
(17) of ref. [39] as 



^(l)SM 



4 top Kh^t 



El{x) 



C7 'tip {f^t) — ~7J^0,l(^) ' ^8 'top {^^t) 



(D8) 
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where x = [m^'^ {^t) / m^Y . To compute x, the MS-renormahzed running top mass is extracted 
from the pole mass according to 



MS/ pole\ r 

rrif. [nil ) = rrif 



4 asjmt 

3 TT 



-, -1 



(D9) 



and is then evolved to the renormalization scale fit according to eq. ( lD4p . 

The new-physics contribution to K can be obtained by simply replacing the top contribu- 
tions in eq. (]D7|) with the corresponding NP contributions. In the MSSM one has to add also 
the squark-gluino contributions to the Wilson coefficients Ci and C2, presented in the appendix 
A. In summary, 



Air 



k=l 



Vt 



(DIO) 



where the "magic numbers" hk are the same as in Kc , and the fk are: 
fk = (0.5784, -0.3921, -0.1429, 0.0476, -0.1275, 0.0317, 0.0078, -0.0031). (Dll) 

In the determination of the electroweak contribution ^ew that appears in eq. flD2p we adopt 
the results of ref. [29] [see in particular eq. (3.10) there], including the available NP contributions 
where appropriate. 

The gluon-bremsstrahlung contribution to the branching ratio is 



(O)eff/ 



TT 



(D12) 



where: 6 = 1 — 2Eo/ml^; the expression for the functions 4>ij{S) are given in eqs. (E.2)-(E.8) 
of ref. [20]; the ratio rrib/ms appearing in the expression of (pssi^) is taken as input. The 
contributions corresponding to3<'j<6or3<j<6 are negligible and we set them to zero. 
We take the SM contributions to the LO effective Wilson coefficients C^'^^^^{nb) from eq. (E.9) 
of ref. [20], and include the NP contributions where appropriate. 

Finally, the non-perturbative contribution to the branching ratio can be approximated as 
(see [20] for the complete list of references) 



N{Eo 



_6_ _ 12 



A, 
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(D13) 



Where the parameter A2 is taken as input, and for the charm quark mass we use {rric). 
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PROGRAM SUMMARY 

Manuscript Title: SusyBSG: a fortran code for BR[i3 Xg 7] in the MSSM with Minimal Flavor 
Violation 

Authors: G. Degrassi, P. Gambino, P. Slavich 

Program Title: SusyBSG 

Journal Reference: 

Catalogue identifier: 

Licensing provisions: None 

Programming language: Fortran 

Operating system: Linux 

Keywords: Supersymmetry, B physics, rare decays 
PACS: 12.60.Jv, 13.20.He 

Classification: 11.6 Phenomenological and Empirical Models and Theories 
Nature of problem: 

Predicting the branching ratio for the decay B ^ Xs^ in the MSSM with Minimal Flavor Violation. 
Solution method: 

We take into account all the available NLO contributions, including the two-loop gluino contributions 
to the Wilson coefficients computed in ref. [1]. The relation between the Wilson coefficients and the 
branching ratio for B ^ Xg'j is computed along the lines of ref. [2] . The SUSY input parameters can 
be read from a spectrum file in the SLHA format. 
Restrictions: 

The results apply only to the case of Minimal Flavor Violation. 
Unusual features: 

Numerical instabilities may arise for special combinations of the input parameters. They can usually 
be avoided by compiling the code in quadruple precision. 
Running time: 

Approximately 10 ms in double precision, 0.5 s in quadruple precision (on a Pentium D 3.40-GHz). 
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